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Abstract. - We focus on the role played by the node degree distribution on the way collective 
phenomena emerge on complex networks. To address this question, we focus analytically on a 
typical model for cooperative behaviour, the Majority Rule, applied to dichotomous networks. 
The latter are composed of two kinds of nodes, each kind i being characterized by a degree 
ki. Dichotomous networks are therefore a simple instance of heterogeneous networks, especially 
adapted in order to reveal the effect of degree heterogeneity. Our main result are that degree 
heterogeneity affects the location of the order-disorder transition and that the system exhibits 
non-equipartition of the average opinion between the two kinds of nodes. This effect is observed 
in the ordered phase and in the disordered phase. 
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Introduction. — It is more and more common, nowa- 
days, to use models and tools from Statistical Physics in 
order to describe the emergence of collective phenomena in 
social systems. Amongst other examples, one may think of 
opinion formation [1,2], rumour or disease spreading [3,4], 
language dynamics [5], etc. In these models, agents are 
located at the node of a graph and are endowed with a 
finite number of available states, e.g. two states - spin up 
and spin down. The links between nodes represent the 
relations (e.g. friendship, co-authorship) between agents. 
Most of the models are based on local attractive interac- 
tions, i.e. agents copy the behaviour of their direct neigh- 
bours, while random changes also take place, thereby mim- 
icking the effect of individual choices. Contrary to classi- 
cal problems of Statistical Physics, though, the underlying 
network structure is not a c?-dimensional regular lattice, 
but a complex network [6], whose node degree (number 
of links per node) and other internal properties may vary 
from one node to another. 

Several works have revealed that a given model may 
exhibit very different (even qualitatively) behaviours de- 
pending on its underlying topology [4]. Important fac- 
tors are for instance the degree heterogeneity [7] , the pres- 
ence of communities [8,9] or the small- world property [10]. 
From a practical point of view, it is therefore important 
to elucidate how these structural properties affect critical 
behaviour if one intends to reproduce the emergence of 



collective phenomena as they take place in realistic sit- 
uations. A typical example would be the propagation of 
rumours in a social network, which is of primordial interest 
for viral marketing issues [11]. From a theoretical point 
of view, it is challenging to develop tools and methods in 
order to determine the influence of the network topology 
on the way the system orders (or not). Such theoretical 
studies have been performed in some specific cases, e.g. 
the Voter Model [7] or the Unanimity Model [12], but a 
coherent and unifying framework is still needed. 

In this Letter, we address this problem by focusing on a 
variant of the majority rule model (MR), that is a typical 
model for consensus formation [1]. Its simple ingredients 
allow a very complete analytical description in the mean- 
field approximation [13]. One should also stress that, con- 
trary to the voter model [14], MR does not conserve av- 
erage magnetization [15]. In the following, we focus on a 
variant of MR that includes random flips (thermal fluc- 
tuations) and study analytically the effect of the degree 
heterogeneity on the phase diagram. To do so, we intro- 
duce dichotomous networks, in which there are two sorts of 
nodes, each sort being characterized by a degree ki or ^2- 
It is shown that the system undergoes a transition from 
a disordered phase to an ordered phase for weak enough 
random effects (~ low temperature) . Our main results are 
that the location of this transition depends on the degree 
heterogeneity 7 = ^2/^1. Moreover, the system is shown 
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to exhibit non-equipartition of the average magnetization 
(each sort of nodes is characterized by a different average 
opinion/spin) when 7^1. This is observed in the ordered 
and in the disordered phase. 

Majority Rule. — The network is composed of A'' 

nodes, each of them endowed with an opinion Oi that can 
be a or (3. At each time step, one of the nodes is ran- 
domly selected. Two processes may take place, i) With 
probability q, the selected node s randomly changes its 
opinion: 

Os a with probabiUty 1/2, 

Og ^ (3 with probabiUty 1/2. (1) 

ii) With probability 1 — q, two neighbouring nodes of s are 
also selected and the three agents in this majority triplet 
all adopt the state of the local majority. The parameter 
q therefore measures the competition between individual 
choices, that have a tendency to randomize the system, 
and neighbouring interactions, that tend to homogenize 
the opinion of agents. In the case = 0, it is well-known 
that the system asymptotically reaches global consensus 
where all nodes share the same opinion [13]. In the other 
limiting case q = 1, the system is purely random and 
the average (over the realizations of the random process) 
number of nodes with opinion a at time t, denoted by At, 
goes to N/2 for t large. 

Homogeneous network. — In this section, we as- 
sume that the network of individuals is highly connected 
and homogeneous, i.e. all the nodes have the same degree. 
In that case, the mean-field rate equation for At reads 

At+i =At + q{^ -a)- 3(1 - q)a{l -3a + 2a^), (2) 

where at = At/ N is the average proportion of nodes with 
opinion a. The term proportional to q accounts for the 
random flips and the last term for local majorities. This 
comes from the fact that the probability for two nodes 
a {(3) and one node /3 (a) to be selected is 3a^(l — a) 
(3a(l — a)^), so that the total contribution to the evolution 
of At is 

W = 3 {a^{l - a) - a(l - af) = -3a(l - 3a -h 2a^). (3) 

Let us stress that Eq. (3) differs from Eq. (2) of [12] by 
a factor 3. In [12], this factor could be absorbed in the 
time scale as it did not play a relevant role. It is straight- 
forward to show that a = 1/2 is always a stationary solu- 
tion of Eq. (2), as expected from symmetry reasons. This 
is obvious after rewriting the evolution equation for the 
quantities A = A - N/2 and S = a - 1/2 

At+i = At+^-{3-5q- 12(1 - q)d^) , (4) 

from which one finds that the symmetric solution a = 1/2 
ceases to be stable when q < 3/5, and that the system 



reaches the following asymmetric solutions in that case 

_ 1 I 3~5q _ 1 / 3-5g 

'^--2-\/l2(T3^' "+-2+VT2(T3^- 

The system therefore undergoes an order-disorder transi- 
tion at q = 3/5. Under this value, a collective opinion 
has emerged due to the im,itation between neighbouring 
nodes. Let us stress that Eqs. (5) respectively converge 
to a_ = and a+ = 1 in the limit g — > 0. 

Degree heterogeneity. — The main goal of this Let- 
ter is to understand how the degree distribution of the 
underlying network affects collective behaviours as those 
of the previous section. Contrary to more phenomeno- 
logical approaches, where the effects of the degree hetero- 
geneity are brought to light by comparing the behaviour 
of a model on several kinds of networks [16] (e.g. Erdos- 
Rcnyi, scale- free, etc), we prefer to address the problem 
from an analytical and more fundamental point of view. 
To do so, we generalize homogeneous networks in the most 
natural way by considering random networks whose nodes 
may be divided into two classes, the nodes in different 
classes being characterized by a different degree, fci or ^2. 
This binary mixture, that we call a dichotommis network 
is particularly suitable in order to reveal the role of degree 
distribution. Indeed, the degree heterogeneity is tunable 
through the parameter 7 = k2/ki. When 7^1, one 
recovers an homogeneous network. Let us also emphasize 
that dichotomous networks differ from usual bipartite net- 
works [17] by the fact that links may exist between nodes 
of the same class. 

The underlying topology is therefore determined by the 
values ki and ^2, and by N^.^_ and N^.^ that are the num- 
ber of nodes of each class. In the following, we will as- 
sume that Nki = Nk2 for the sake of simplicity. A more 
complete analysis for the full range of parameters will be 
considered elsewhere. We are interested in Ai-t and A2;t 
that are the average number of nodes 1 and 2 with opin- 
ion a. By construction, each node is selected with the 
same probability during one time step, but those with a 
higher degree have a larger probability to be among the 
neighbours of the selected node, i.e. to be in the majority 
triplet. This effect will have to be taken into account in 
order to generalize Eq. (2). 

By construction, the probability that the selected node 
has a degree ki is 1/2, but the probability that the neigh- 
bour of this selected node has a degree kj is kj/{k\ + fe) 
(one assumes that there are no correlations between the 
degrees of neighbouring nodes). Consequently, the proba- 
bility that the selected node has a degree ki and that both 
of its selected neighbours have a degree k2 is 

2(/ci + /c2)2 2(1 + 7)2- ^ > 

Similarly, the probability that the selected node has a de- 
gree /ci, that one of its neighbours has a degree fci and 
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that the other neighbour has a degree ^2 is 

fcifc2 _ 7 
{h+hr ~ (1 + 7)'' 



(7) 



while the probability that all three nodes have a degree ki 
is 



1 



2(fci + fc2)2 2(1 + 7)2- 



(8) 



The sum of these three probabilities is normalized and the 
probabilities of events when the selected node has a degree 
k2 are found in the same way. Putting all contributions 
together, one finds the probabilities P(x,y) that x nodes 1 
and y nodes 2 belong to the majority triplet 



^"(3,0) = 



^(1,2) = 



1 



2(1+7)^' 
27 + 72 

2(1+7)2' 



_ 1 + 27 

7' 



(0,3) 



2(1 + 7)2 



(9) 



where the normalization '^^y P{x,y) = 1 is verified. In or- 
der to derive the coupled equations generalizing Eq. (2) 
for Ai-t, one needs to evaluate the possible majority pro- 
cesses taking place when a triplet {x, y) is selected. Let us 
focus on the case (2, 1) as an example. In that case, the 
number of nodes A\, A2 will change due to the contribu- 
tions 



1,(2,1) 



[20102(1 - oi) - 2ai(l - a2)(l - ai)] 



2,(2,1) 



[0^(1 - 02) - 02(1 - Oi) 



(10) 



where is the proportion of nodes with opinion a in the 
class i. The first line accounts for cases when one node 1 
and one node 2 have the same opinion but disagree with a 
node 1, while the second line accounts for cases when the 2 
nodes 1 have the same opinion but disagree with the node 
2. The other situations {x,y) arc treated similarly [19]. 
Putting all the contributions together, one arrives at the 
set of non-linear equations 



Ai 



■,t+i 



. q qai 

- M;t = T TT 



[^a\bi - aibl) 



A 



2;t+l 



l2;t 



4 2 ' 2(1+7)2 

+ 2(l + 27)(a2ai&i-ai&2?>i) 

+ {2^ + j^)(albi-a,bl)] 

q qa2 , (1 - (?) 2/ 2, .2n 
= i" — + ^(TT^[37(a262-a2M 

+ 2(27 + 72)(aia2&2 - 026162) 

+ {l + 2j){a%-a2bl)], (11) 



where bi is the proportion of nodes with opinion (3 in the 
class i (bi = 1 — ai). One verifies by summing the equa- 
tions for Ai and A2 that Eq. (2) is recovered in the limit 
7=1 (with a = (ai + a2)/2), as expected due to the 
indistinguishability of the nodes in that case. 

It is easy to show that ai =02 = 1/2 is always a station- 
ary solution of the above set of equations. The stability 
analysis is performed by looking at deviations to this so- 
lution ai = 1/2 + 5i, 02 = 1/2 + ^2 and keeping only 
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Fig. 1: Location of the order-disorder transition qc as a function 
of 7. The x-scale is logarithmic in order to make the symmetry 
7 <-> 7~^ clear. 

linear corrections. In the continuous time limit and after 
rc-scaling the time units, the evolution equations for these 
deviations read 

„ . -3 + 5g + 27(l + g)+72(l+g) 

4(1 + 7)2 
(l + 47 + 72)(l-g) ^ 

+ WT^' 

o . l + g + 27(l + g)+72(-3 + 5g) , 
Otd2 = T7:r~. ^ "2 



4(1 + 7)2 
(l + 47 + 72)(l-g) ^ 

+ WTW' 



(12) 



The stability analysis requires to look at the eigenvalues of 
the above evolution matrix and to find the critical value of 
q at which the real part of one eigenvalue becomes positive. 
Lengthy calculations lead to 



9c (7) = 



- I67 + 872 + 16VK 



24 + I67 + 2472 + 16VK 



(13) 



where K = 2 + 8-f + I672 + 87^ + 27^^. In the limiting 
case 7-^1, one recovers the known result qc = 3/5. 
It is also possible to verify (see Fig. 1) that this re- 
lation is symmetric under the changes 7 <-> 'y~^, i.e. 
under an exchange of nodes 1 and 2. Moreover, the 
maximum value is obtained for 7 — > and 7 — > 00, 
9c(0) = gc(oo) = (1 + 2\/2)/(3 + 2\/2). Our first con- 
clusion is therefore that the location of the order-disorder 
transition depends in a non-trivial way on the degree het- 
erogeneity 7, even though these deviations remain of small 
magnitude (maximum 10%). Moreover, the location of the 
transition is shifted to higher values of q when the system 
is more hereogencous and qd^) < qda). 

It is important to remind that node degree heterogene- 
ity leads to a heterogeneity of the probability to belong to 
a majority triplet (sec Eqs. (9)) and therefore to a het- 
erogeneity of the frequency of interaction. In contrast, the 
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frequency of random flips remains the same for the two 
types of nodes. Howcv(^r. as emphasized above, the; order- 
disorder transition comes from a competition between the 
frequency of random flips ~ q and the frequency of major- 
ity triplets ^ (1 — (j). It is therefore understandable that 
the heterogeneity of the system might lead to a shift of Qc- 
In order to highlight this effet in detail, it is interesting 
to consider the limiting case 7 = 0, where the network is 
so heterogeneous that any randomly chosen link arrives at 
a node of type 1. Consequently, a node 2 has a vanish- 
ing probability to be chosen as a neighbour of the selected 
node and it may belong to a majority triplet only if it is 
the selected node itself. When 7 = 0, the equations of 
evolution for the quantities Ai read 



t+1 



Av,t 



q qai ^ (1 - g) 2, .2n 
4 - ^ + — 2 — t^'^i^ ~ ^' 



i2;t+l 



-I- 2(a2ai6i - 0162^1)] 
= 9 _ 902 (1 - q) 

4 2 2 



[(af 62 - a2bl) 



and it is straightforward, but lengthy, to show that the 
stationary solutions are either ai =02 = 1/2 or ai = 
l/2±5i,ai = 1/2 ±(52 with 



Si 
S2 



^2(g2 - 1) + (1 - 9) v/25^^22^T? 



Si 



12(1 -g) 



(14) 



It is also possible to show that the ordered solutions exist 
only for q < qdO), where qdO) is given by (13). It is 
instructive to look at the linearized dynamics of Si around 
the disordered solution 

dtSi = (3 - 5q)Si + 2(1 - q)S2 

dtS2 = 2{1 - q)Si - {1 + q)S2, (15) 

whose matrix of evolution Lij is obviously asymmetric un- 
der the change 1^2, contrary to the case 7 = 1, where 
the linearized equations read 



dtSi = -4qSi + 6(1 - q)S2 
dtS2 = 6(1 - q)Si - AqS2- 



(16) 



These two sets of equations not only differ by their sym- 
metry, but also by the relative importance of the diag- 
onal terms La as compared to the cross terms Lij i ^ j 
(such terms measure the coupling between the two kinds of 
nodes). This difference makes possible a shift of the value 
of qc where the determinant of Lij vanishes, i.e. where 
one eigenvalue vanishes. 

Let us now return to our study of MR for general val- 
ues of 7. In order to elucidate the behaviour of ai and 
a2 below qc, we have performed numerical integration of 
Eqs. (11). It appears (see Fig. 2) that ai and 02 reach 
different asymptotic values ai^oo 7^ 02,00 and that the class 
of nodes with the higher degree exhibit larger deviations 



Ys=l0,5i 
Ys=lO, 82 
Numerical integration 




1.6 
1.5 
1.4 
1.3 
1.2 
1.1 





Yr=10 □ 


Numerical Integration 


Theoretical prediction 




□ 








□ \ " 




□ \ 




□ 




□ - 




□ 




□ 




□ 



0.2 



0.4 



0.6 



0.8 



Fig. 2: Bifurcation diagram of 5i{q) (upper figure) and R{q) 
(lower figure) for 7 = 3.162 (JV = 10* nodes, pi = 0.005, 
7s = 10, see section Simulation results). Solid lines are 
the asymptotic solutions obtained by integrating numerically 
Eqs. (11). Simulation details arc found in the main text. In 
the lower figure, the dashed line correspond to the theoretical 
predictions Eq. (17) and Eq. (19). The vertical line indicates 
the value of q above which S is measured during the relaxation 
to Si = 0, and not anymore in the asymptotic state 6i ^ 0. 



to 1/2 than the other class. This may be understood 

by the fact that nodes with a higher degree are more of- 
ten selected in majority triplets, thereby triggering their 
tendency to reach consensus. In order to evaluate this 
non-equipartition of the average opinion between the two 
species of nodes, we introduce the ratio R = S2/S1. One 
observes (not shown in Fig. 2) that R = 1 when 7=1 
and that increasing values of 7 lead to increasing values 
of R at fixed q. 

We tackle this problem from an analytical point of view 
by focusing on the limit g — > and looking for small de- 
viations to global consensus ai = + Ci. By inserting this 
development into Eqs. (11) and keeping linear terms, it 
is straightforward to show that the asymptotic values of 
ai are ai = 2(^1) 1 ^2 = 2(1+57) ^' from which one 
shows that 
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Fig. 3: Typical realization of a dichotomous network composed 
of = 50 nodes, with pi — 0.05 and 7s = 5. In the lower left 
corner, we plot the histogram of the node degrees that show two 
peaks. Each peak is associated to one class of nodes. The net- 
work is visualized thanks to the visone graphical tools [21] with 
a radial layout, i.e. the higher the degree of a node, the more 
central is its position in the graph. The hidden variable het- 
erogeneity obviously leads to degree heterogeneity, while pre- 
serving its random structure. 
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Fig. 4: Histogram of the node degrees for a dichotomous net- 
work composed of A'^ = 10* nodes, with pi — 0.005 and 
7s = 10. The vertical lines indicate the average degrees 
ki ^ 104 and « 329 obtained from Eqs. (20). By con- 
struction, the surface under each sub-distribution is equal to 
N/2 = 5 X 10^ 



(18) 



i?w H-4 



7 



(5 + 7)(l + 57)- 



(17) whose stationary solution is found to be 



This solution is in perfect agreement with the numeri- 
cal integration of Eqs. (11) in the limit of small q, and 
is asymmetric, i.e. R{j~^) = R{'y)~-^ (which leads to 
R{'y~^) — 1 = — (i?(7) — 1) in the linear approximation) 
under the change 7 ^ 7^^, as expected. One can also 
show that Eq. (17) with 7 = is recovered from the exact 
solution (14), by dividing 62 by di and keeping only first 
corrections in q. 

Before going further, we would like to focus on a prop- 
erty that may appear intriguing at first sight. Namely, 
if one calculates the time evolution of R from Eqs. (11) 
above the critical point qc, one finds that this quantity 
asymptotically converges to a constant value different from 
zero Roo 7^ 0. In contrast, we have shown above that 
the deviations Si and S2 go to zero for these values of q 
(ai = 1/2 is a stable solution). Altogether, this suggests 
that the relaxation to the asymptotic state behaves like 
^1 — Cie^'^i*, 62 = 6*26^'^^*, where the relaxation rates 
are the same Ai — X2 and where C2 = RooCi. This be- 
haviour reminds the non-equipartition of energy taking 
place in inelastic gases, for which it is well-known that 
different components of a mixture may exhibit different 
partial temperatures [18]. This analogy is clear after not- 
ing the similarities between the linearized equations (11) 
and the equations for the partial temperatures Ti [19] ob- 
tained from Inelastic Maxwell Models [20]. 

By using the same methods as in [19], one writes a closed 
equation for R 



dtR 



2(1 + 7)' 



•[(I-I-47 + 72) -f 2(72 - 1)R 



R 



(7^ - 1) + ^/2^1 + A-f + 872 + 473 + 74 
(1 -f 47-1-72) 



(19) 



This relation shows the same asymmetry as Eq. (17). One 
also verifies that R goes to 1 when 7^1 and that R goes 
to the finite values (1-1- \/2)~^ and (1 -I- \/2) when 7^0 
and 7 ^ 00 respectively. 

Simulation results. — By construction, the random 
steps of MR are easy to implement in a computer sim- 
ulation. The only difficulty is to build the underlying 
dichotomous network, i.e. a network in which there are 
two characteristic node degrees but without internal cor- 
relations. We tackle this problem by applying a method 
reminiscent of networks with hidden variables [22-24] and 
Superstatistics [25-27]. Namely, we first prepare N nodes 
and assign them hidden variables hi < 1, half of the nodes 
receiving the hidden variable hi = pi and the other half 
the hidden variable hi = P2- Next, each pair of nodes («, j) 
receives a link with a probability equal to yjhihj < 1 (see 
Fig. 3). Let us introduce 7s = P2/P1 that measures the 
hidden variable heterogeneity. From the above definitions, 
one finds that the degree distributions of the nodes 1 and 
2 are peaked around 

N , 



ki = '—{pi + VP1P2) 
k2 = y(p2 + VPiP2)- 



(20) 



The resulting network (see Fig. 4) is therefore a good 
candidate for dichotomous networks, where the effective 
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7 = fc2/A;i is related to 75 through the relation 



^ P2 + _ 75' + 

Pi + VpTP2 1 + x/Ts 



(21) 



Let us stress that the peaks of the degree distribution are 
not exactly delta functions, but their width is vanishingly 
small for large enough values of ki and k2- 

The computer simulations presented in this Letter are 
performed with TV = 10^ nodes, pi — 0.005 and 75 ~ 10, 
but other parameters have also been considered and lead 
to the same conclusions. The corresponding degree het- 
erogeneity is therefore 7 w 3.162. The stationary state of 
MR is measured by starting the simulation from the initial 
condition Ui = 1/2 and letting the system relax toward its 
asymptotic state during 100 time steps/node. Afterwards, 
one measures \Si \ and \S2\ during 500 time steps/node, av- 
erage over the time and over 100 realizations of the random 
process. This method gives an excellent agreement (see 
Fig. 2) with the numerical integration of Eqs. (11) below 
the critical value Qc, but it is not applicable for measuring 
R above this critical value. This is due to the fact that 
non-equipartition of opinion takes place during the relax- 
ation to the asymptotic solution Oj = 1/2 when q > Qc 
and that this relaxation is indiscernible from the finite 
size fluctuations around = 1/2 in the long time limit. 
Consequently, we use another strategy in order to measure 
R when q > qc- Namely, the simulation is started from 
tti = 0.7, but R is now evaluated during the relaxation, i.e. 
after 50 time steps/node, and averaged over 5000 realiza- 
tions of the random process. These simulation results arc 
in very good agreement with Eq. (19) (at least for q < 0.8, 
above which fluctuations are very large) and confirm that 
R docs not depend on q in the disordered phase. Finally, 
let us stress that the parameters used for the simulations 
imply that the degree of the nodes is very large (see Fig. 
4). This choice ensures that the diameter of the network 
is very small and that the mean-field approximation used 
in order to derive Eqs. (11) is justifiable. It is therefore 
expected that simulations performed with lower values of 
pi exhibit discrepancies from the theoretical predictions. 

Short discussion. — In this last section, we would 
like to discuss the implications of our theoretical predic- 
tions. The non-equipartition of opinion observed in MR 
is a striking feature that implies that the state of a so- 
cial agent strongly depends on its connectivity. One ex- 
pects that this effect should also take place in other models 
for opinion formation, at least if the frequency of interac- 
tion between agents is heterogeneous. Non-equipartition 
of opinion should be searched in empirical data, in the 
Blogosphere for instance, by averaging the state of agents 
having the same connectivity and looking for a relation 
between this average value and the connectivity. Such em- 
pirical study could also useful in order to test the validity 
of MR as a relevant model for interacting social agents. 

In this paper, we have also shown that the location of 
the order-disorder transition depends on the degree het- 



erogeneity. This shift might have striking consequences in 
realistic; situations as it implies that a change of the un- 
derlying topology might lead to a transition and to the 
emergence of a new (ordered or disordered) phase in the 
network. Similar shifts should also be searched in other 
models for cooperative behaviour, e.g. Ising models, lan- 
guage dynamics, etc, when they are applied to heteroge- 
neous networks and, in particular, to dichotomous net- 
works. 
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